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1. Introduction 



In the past few years many new analytical results have emerged from Matrix Models (MM) 
for the Dirac operator spectrum with non- vanishing chemical potential /I / 0. There are 3 different 
possible chiral symmetry breaking patterns [JXJ] and in the MM picture they differ crucially in the 
way the eigenvalues are depleted from the imaginary axis through ji. ^ [§]. Today we have 
detailed predictions for microscopic Dirac spectra in 2 of these symmetry classes, both quenched 
and unquenched: QCD and the adjoined representation, where the latter is replaced by SU (2) 
colour in the fundamental when using staggered fermions as we will do. After first approximate 
results for quenched QCD [^] the exact quenched density was derived in [Q] and related to chiral 
Perturbation Theory in the epsilon regime (e#PT). Unquenched partition functions were computed 
in [||] and related to £#PT, and finally fully unquenched Dirac spectra for QCD became available 
[§,0]. Results for the adjoint (or SU (2) staggered) class followed very recently ^ and we confront 
its unquenched predictions including dependence on /J. and quark mass m with Lattice data here. 

Not all the above MM results have been compared to the Lattice so far, precisely due to the 
sign problem in unquenched QCD. Up to now only quenched simulations at \i / using staggered 
fermions have been successfully described: for QCD []|] (see also [[K]] in these proceedings) and 
for SU (2) [pT|], which we will extend here. Previous comparisons fli2"| ] were lacking analytic pre- 



dictions at the time, they could only be done in the bulk of the spectrum [13] for MM without 
chiral symmetry. Very recently it has been shown at /x = that the previous topology-blindness of 
staggered fermions can be cured by improvement, as reviewed in [14]. 

In the last years different ways of attacking the sign problem in QCD were developed: multi- 
parameter reweighting, Taylor-expansion and imaginary jj, (see Jl5| ] for a review and references). 
However, none of these have been applied so far to the region of £#PT where a comparison to MM 
is expected to hold. We purse a different avenue here, choosing an SU (2) gauge theory without sign 
problem where dynamical simulations can be performed in a standard way []T(J. This permits us to 
extend previous MM comparisons Jl7| , |l8| ] to jj. ^ 0. It is of principle interest to test the validity of 
all MM predictions for complex Dirac spectra on the Lattice including dynamical fermions. 



2. Predictions from Matrix Models with /i ^ 

In this section we briefly introduce the relevant MM used and give its results for the spectral 
density. For more details and references we refer to [Q]. The MM partition function is given by 

&$"'\{m f }\pL)-= Jdd>dVexp[-NTr(&<l> + W)] fl det (^t7^t ^ij' {2A) 

where 1 is the quaternion unity element. The two rectangular (N + v) xN matrices, <1> and *P, 
contain quaternion real elements and replace the off-diagonal blocks of the Dirac operator, averaged 
with a Gaussian matrix weight instead of the gauge action. Here we model the jiiyo part with a 
random matrix of the same symmetry as the kinetic part, assuming it is non-diagonal in matrix 
space (as in [||] for QCD). If universality holds this choice of basis should not matter, compared 
to *P replaced by unity as [Q]. This assertion is true in the QCD symmetry class, see [Q, Q] vs. 
[|], [7|]. The size of the Dirac matrix 2Af ~ V relates to the volume, and we have chosen rectangular 
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matrices to be in the fixed sector of v > zero eigenvalues or topological charge. Because of using 
unimproved staggered fermions we set V = in our comparison to the data later. 

After transforming the linear combinations £<E>W + jU^F^ to be triangular we can change to 
complex eigenvalues Zk=i,...,N (rotated to lie on the real axis for pt = 0). All their spectral correlation 
functions can be computed using skew orthogonal Laguerre polynomials in the complex plane [g]. 
In the large-Af limit the complex eigenvalues, masses and chemical potential have to be rescaled 



y/l Nz = £, , V2Nm = rj, and 



lim 2N/A* 



or 



(2.2) 



This limit is called weakly non-hermitian as /I 2 is rescaled with the volume. The same scaling is 
found in the static part of the chiral Lagrangian [ 19] dominating at £#PT. While the latter contains 
two parameters, the chiral condensate £ and the decay constant F n , the MM contains only one, the 
variance fixed to be £ = yl here. Thus the constant multiplying a is not know and we are left with 
1 free fit parameter, in contrast to the parameter-free MM prediction at pt = 1 . 

The unquenched microscopic spectral density, normalised to 3m(^)5'(3/M(^)) for a — > 0, 
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splits into the quenched part and a correction term Ap^ ak (^ ; rj) depending on the mass rj 
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We only give the result for Nf = 2 staggered flavours of equal mass here, see [g] for more flavours 
and higher correlation functions, as well as [17] for matching MM and staggered flavours. Fig. 
[T] shows this density in the complex plane at weak non-Hermiticity, as we shall compare to data 
below. The eigenvalues are repelled from the real axis for ji / 0, being a distinct feature of this 
symmetry class (compare to [^] for QCD). We have observed this repulsion in the data for pi as 
small as 10~ 6 . The individual eigenvalues located at the maxima previously are now split into a 
double peak in the complex plane. 

Increasing the mass 17 — > 00 moves the density to the left, bringing it back to the quenched 
expression, eq. (2.3) at &p^ ak (E,;r\) = 0. Decreasing r\ — » pushes the eigenvalues further away 
from the origin, approaching the quenched density at v = 2 approximately. Increasing a rapidly 
washes out the oscillations and leads to the formation of a plateau. The limit a — > °° takes us to the 
MM at strong non-Hermiticity, with pi unsealed (see [Jsj] ) - A comparison to quenched Lattice data 



in this regime was given previously in ply, our unquenched results will be reported elsewhere. 
'This observation was missed in |p[[llp. 
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Figure 1: The Nf = 2 flavour spectral density at a = 0.012 with masses rj = 8.74 and v = (left) and a cut 
through the density normalised to unity for different mass values compared to quenched (right). 



3. Lattice data with dynamical fermions at /i ^ 

Our data were generated for gauge group SU (2) with coupling j8 = 4/g 2 = 1.3 and Nf = 2 



staggered flavours, using the code of Jig]. In this setup the fermion determinant remains real (as it 
does in the MM, see [JsJ] ) and standard Monte Carlo applies. We have studied two different Volumes 
6 4 and 8 4 for various values of /I = 10~ 6 — 0.4 and values of the quark masses ma = 0.025 — 20. 
Because the eigenvalues lie in the complex plane of the order of 5- 10k configurations are needed. 
In order to have a window where a MM description applies for these small lattices we have to go 



to relatively strong coupling (see [10] for a discussion of this issue in QCD). 



To compare with the prediction ( |2.3| ) we have taken cuts through the density: along the maxima 
parallel to the real axis, Fig. [I] right, and perpendicular to that over the first maximum pair. The 
effect of dynamical fermions is most clearly seen in the shift in the first cut where we choose the 
values rj = 14 (blue) and T] = 33 (pink) to be used in our data below, compared to quenched (black). 
Being very costly we did not to go to smaller t]. 

The parameters rj and a were obtained as follows from our data with input ma and }i. First, 
we determined the rescaling of the masses and eigenvalues by measuring the mean level-spacing 
d ~ l/p(0), using the Banks-Casher relation from jU = 0: 7ip(0) = LV where p(0) is the mean 
spectral density. Due to jU <C 1 the geometric distance between eigenvalues agrees within errors 
with the distance obtained by a projection onto the real axis. This provides us with the rescaling of 
the eigenvalues and masses, 

zan/d =<§ and man/d = rj . (3.1) 

At the same time the spacing d contains the volume factor for the rescaling of ji 2 : a 2 = C\l 2 njd. 
The constant C ~ F 2 /L of order 1 is obtained by fitting the data to the cut parallel to the imaginary 
axis on the first maximum. Since the surface under this curve is finite we can fit to its integral 
function, being independent of the choice of histogram widths in imaginary direction. This also 
fixes the normalisation. For illustration purposes we show the cut and not its integral. We can now 



test the scaling hypothesis of ji 2 eq. (|L2|). For this purpose we have kept ji l V fixed for the two 



2i 



volumes V = 6 4 and 8 4 : jl = 1 • 10 3 and jX = 5.625 • 10 4 , respectively. Since the level spacing 
d also depends on the mass we have kept ma fixed, leading automatically to different rj -values for 
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different volumes 2 . The data (histograms) in Fig. § right confirm this scaling very well: the fitted 
a = 0.013 from V = 6 4 (blue curve) describes the V = 8 4 data as well (green). For V = 8 4 we also 
display the a = 0.012 (blue) obtained from an independent fit, they agree within 5%. In the cuts 
in Fig. ^| left these small variations in a cannot be seen. If we were to compare to the quenched 
density a different fit value for a ~ 0.0185 (6 4 ) and a ~ 0.0175 ( 8 4 ) would be obtained instead, 
describing the right curves equally well. However, in the left plots the quenched MM curve (grey) 
deviates from the unquenched (blue) one. In the upper V = 6 4 plot the discrepancy from tj = 33 to 
quenched can be clearly seen in the data as we capture up to the 4th maximum, whereas in V = 8 4 
keeping ma fixed implies rj = 102, taking us back to almost quenched. 

In order to see the difference from quenched at smaller masses we compare to data correspond- 
ing to a rescaled mass tj = 14 (blue curve) in Fig. || below. Although we can only resolve well the 
first 2-3 maxima in the left picture, the mismatch with the quenched curve (grey) is evident. 

To summarise we have shown that the MM correctly predicts complex SU (2) Lattice data with 
Nf = 2 dynamical staggered fermions at pL ^ 0, describing the effect of small quark masses. We 
have also confirmed the scaling of jl 2 with the volume at weak non-Hermiticity from unquenched 
Lattice data. 
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2 For the comparison in the different tj ^> 1 were both close enough to quenched. 
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Figure 3: Data (histograms) vs. MM (blue curve=unquenched, grey=quenched) for V = 6 4 ma = 0.035 and 
ju = 1 • 10~ 3 => 77 = 14 and a = 0.012. 
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